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ABSTRACT 

We calculate the stellar mass-metallicity relation at five epochs ranging to z ~ 2.3. We quantify 
evolution in the shape of the mass-metallicity relation as a function of redshift; the mass-metallicity 
relation flattens at late times. There is an empirical upper limit to the gas-phase oxygen abundance 
in star-forming galaxies that is independent of redshift. From examination of the mass-metallicity 
relation and its observed scatter we show that the flattening at late times is a consequence of evolution 
in the stellar mass where galaxies enrich to this empirical upper metallicity limit; there is also evolution 
in the fraction of galaxies at a fixed stellar mass that enrich to this limit. The stellar mass where 
metallicities begin to saturate is ~ 0.7 dex smaller in the local universe than it is at z ~ 0.8. 
Subject headings: galaxies: abundances — galaxies: ISM — galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

Gas dynamics and star formation regulate the gas- 
phase oxygen abundance (metallicity) of star-forming 
galaxies. Oxygen in the universe forms in the late stage 
evolution of massive stars. Oxygen dispersed into the in- 
terstellar medium (ISM) by supernovae and stellar winds 
increases the metallicity of galaxies as they increase their 
stellar mass. Observations suggest that in most star- 
forming galaxies at z < 2, the increase in stellar mass 
is dominated by cosmological inflows of gas from the in- 
tergalactic medium (Noeske et al. 2007; Whitaker et al. 
2012, among others). At the same time, outflows of gas 
are ubiquitously observed in star-forming galaxies out to 
z ~ 3 (Weiner et al. 2009; Chen et al. 2010; Steidel et al. 
2010). Because metallicity is established by the interplay 
of gas flows and star-formation, observations of the chem- 
ical evolution of galaxies provides important constraints 
for these physical processes in models of galaxy evolution 
(e.g. Dave et al. 2011; Zahid et al. 2012b; M0ller et al. 
2013). 

The mass-metallicity (MZ) relation is a measure of 
the average gas-phase oxygen abundance as a function 
of stellar mass. Lequeux et al. (1979) first observed the 
relation between stellar mass and metallicity in nearby 
irregular and blue compact galaxies. From the analysis 
of ~ 50, 000 star- forming galaxies in the Sloan Digital 
Sky Survey, Tremonti et al. (2004) establish the MZ re- 
lation in the local universe. Subsequent efforts show that 
the MZ relation extends to low stellar masses (~ 10 7 M Q ; 
Lee et al. 2006; Zahid et al. 2012a; Berg et al. 2012) and 
out to at least z ~ 3 (Savaglio et al. 2005; Erb et al. 2006; 
Maiolino et al. 2008; Mannucci et al. 2009; Zahid et al. 
2011; Foster et al. 2012; Yabe et al. 2012; Perez-Montero 
et al. 2012; Yuan et al. 2013, among others), with indica- 
tions that it may extend to even higher redshifts (Laskar 
et al. 2011; M0ller et al. 2013). The metallicities of galax- 
ies increase with stellar mass and the metallicities at all 
stellar masses decrease with redshift. Because of system- 
atic uncertainties and a lack of sufficiently large samples 
of galaxies with well determined metallicities at interme- 



diate and high redshifts, evolution in the shape of the 
MZ relation is not well known. 

A robust determination of the evolution of the MZ rela- 
tion requires consistent measurements of the metallicities 
and stellar masses of galaxies. Here we bring together 
a homogeneously analyzed sample of galaxies with well 
determined MZ relations to investigate the evolution of 
the MZ relation to z ~ 2.3. In Section 2 we describe 
the data and the methods applied to determine metal- 
licities and stellar masses. In Section 3 we present the 
MZ relation and its scatter and quantify its evolution. 
We discuss our results in Section 4 and conclude with a 
summary in Section 5. We assume the standard cosmol- 
ogy (H ,n m ,n A ) = (70 km s" 1 Mpc" 1 , 0.3, 0.7) and a 
Chabrier (2003) IMF. 

2. DATA AND METHODS 
2.1. The Data 

We derive the MZ relation at z = 0.08,0.29 and 
0.78 using data from the Sloan Digital Sky Survey DR7 
(SDSS; Abazajian et al. 2009), the Smithsonian Hec- 
tospec Lensing Survey (SHELS; Geller et al. 2005, Hwang 
et al., in prep) and the Deep Extragalactic Evolutionary 
Probe 2 Survey DR3 (DEEP2; Davis et al. 2003), re- 
spectively. The MZ relation at z = 1.41 and 2.26 are 
determined by Yabe et al. (2012) and Erb et al. (2006), 
respectively. We determine the stellar masses of galaxies 
in their sample using our method (see Section 2.2) and 
adopt their metallicity estimates but convert them for 
consistency with our measurements of metallicity in the 
z < 1 samples (see Section 2.3). 

Our primary selection criteria for galaxies in the z < 1 
samples are the signal-to-noise (S/N) ratios of strong 
emission lines. Foster et al. (2012) show that the de- 
rived MZ relation does not vary significantly with S/N 
cuts on most of the strong emission lines used in this 
study, though S/N cuts on the [OIIIJA5007 line can lead 
to a significant bias. The requisite lines for determin- 
ing metallicity in the z < 1 samples are [OII]A3727, Hf3 
and [OIII]A5007 (see Section 2.3). We require a S/N 
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> 3 in the [OIIJA3727 and H/3 emission lines. Addition- 
ally, in the SDSS and SHELS sample we use the Ha and 
[NII]A6584 emission lines in selecting star- forming galax- 
ies and apply the same S/N criteria. 

The SDSS spectroscopic sample consists of ~ 900,000 
galaxies primarily in the redshift range of < z < 
0.3. We adopt emission line fluxes measured by the 
MPA/JHU 1 group and determine the stellar masses from 
the ugriz-hand photometry. AGN are removed from 
the sample using the [OIIIJ/H/3 vs. [Nil] /Ha diagram 
(i.e. the BPT method, Baldwin et al. 1981; Kauffmann 
et al. 2003; Kewley et al. 2006). We require an aper- 
ture covering fraction > 20%(Kewley et al. 2005) and 
select galaxies with z < 0.12. Our final sample consists 
of ~ 51,000 star- forming galaxies in the redshift range 
of 0.02 < z < 0.12 

The SHELS survey consists of ~ 25,000 galaxies in 
the Fl (Hwang et al., in prep) and F2 (Geller et al. 
2005; Hwang et al. 2012; Geller et al. 2012) regions of 
the Deep Lens Survey (Wittman et al. 2002) spanning 
the redshift range of < z < 0.8. We determine stel- 
lar masses from the ugriz-b&nd SDSS photometry and 
remove AGN using the BPT method. Our sample con- 
sists of 3,577 star-forming galaxies in the redshift range 
0.2 < z < 0.38. 

The DEEP2 2 survey (Davis et al. 2003) consists of 
~ 50,000 galaxies spanning the redshift range of 0.7 < 
z < 1.4. We determine stellar masses from BRI- band 
photometry with additional X s -band photometry for half 
of the sample. This sample differs slightly from the sam- 
ple we analyzed in Zahid et al. (2011). For consistency 
with the SDSS and SHELS analysis, we use stellar pop- 
ulation synthesis models to remove Balmer absorption. 
Additionally, we limit AGN contamination by removing 
17 x-ray galaxies in the sample (Goulding et al. 2012). 
Our final sample consists of 1,254 star-forming galaxies 
in the redshift range 0.75 < z < 0.82. 

Yabe et al. (2012) and Erb et al. (2006) measure the 
MZ relation at z ~ 1.4 and z ~ 2.3, respectively. Yabe 
et al. (2012) determine the MZ relation from the stacked 
spectra of 71 objects with significant detections of Ha 
and Erb et al. (2006) determine the MZ relation from the 
stacked spectra of 87 UV-selected star-forming galaxies. 
The MZ relations determined by Erb et al. (2006) and 
Yabe et al. (2012) are subject to some systematic uncer- 
tainties (see Section 2.3) that we are currently unable to 
quantify. We include these data to demonstrate that the 
best data currently available at high redshift are qualita- 
tively consistent with the main conclusions of this study 
based solely on the z < 1 samples. 

2.2. Stellar Mass Determination 

We measure stellar masses for the five samples us- 
ing the Le Phare 3 code developed by Arnouts & Ilbert. 
We determine the mass-to-light ratio by fitting the SED 
with stellar population synthesis models and we scale 
the luminosity by the mass-to-light ratio to yield a stel- 
lar mass estimate (see Bell et al. 2003). We use the 
stellar templates of Bruzual & Chariot (2003) and a 

1 http:/ /www.mpa-garching.mpg.dc/SDSS/DR7/ 

2 http://deep.ps.uci.edu/DR3/ 

3 http : //www.cf ht.hawaii.edu/-arnouts/LEPH ARE/cfht_ 
lephare/ 'lephare.html 



Chabricr (2003) IMF. The models have two metallici- 
ties and seven exponentially decreasing star formation 
models (SFR cx e -*/ T ) with t = 0.1,0.3,1,2,3,5,10,15 
and 30 Gyrs. We apply the extinction law of Calzetti 
et al. (2000) allowing E(B-V) from to 0.6 and a stellar 
population age range from to 13 Gyrs. For each of our 
samples, the age never exceeds the age of the universe at 
the median redshift of the sample. 

We adopt the median of the mass distribution as our es- 
timate of the stellar mass. We compare the stellar masses 
we determine for our SDSS sample with those derived by 
the MPA/JHU group. We find a 0.17 dex dispersion after 
correcting for a systematic offset between the two meth- 
ods. This dispersion is consistent with the observational 
uncertainty in the photometry. 

2.3. Metallicity Determination 

Ratios of the strength of collisionally excited emis- 
sion lines to recombination lines are both theoretically 
and empirically calibrated for determining metallicity. 
A long-standing astrophysical problem is that the abso- 
lute metallicity determined from theoretical calibrations 
is typically ~ 0.3 dex higher than metallicities deter- 
mined using empirical calibration methods. However, we 
emphasize that our analysis requires a robust relative es- 
timate of the metallicity, which the various metallicity 
diagnostics deliver (Kewley & Ellison 2008). Kewley & 
Ellison (2008) provide a table of formulae allowing for 
conversion of metallicities into various baseline methods. 

We determine metallicities from the i?23 line ratio cali- 
brated by Kobulnicky & Kewley (2004, hereafter KK04). 
A major advantage of this method is that it explicitly 
solves and corrects for the ionization parameter. The 
relevant ratios of measured emission line intensities are 



i?23 



[OILU3727 + [OIII]A4959, 5007 



and 



032 = 



H/3 

[OIII]A4959,5007 
[OII]A3727 ' 



(1) 
(2) 



We assume that the ratio of the fluxes of [OIIIJA5007 
to [OIILU4959 is 3 (Osterbrock 1989) and adopt a 
value of 1.33 times the [OIIIJA5007 when summing the 
[OIII]A5007 and [OIII]A4959 line strengths. The i?23 
method is sensitive to the ionization state of the gas and 
the 032 ratio is used to correct for variations. We ap- 
ply this method to the SDSS, SHELS and DEEP2 data. 
The galaxies in this study all lie on the upper metallicity 
branch (see Zahid et al. 2011). 

At higher redshifts the requisite lines for determination 
of the i?23 ratio are rarely observed. Pettini & Pagel 
(2004, hereafter PP04) calibrate 



N2 = ([NII]A6584/Ha) 



(3) 



to yield the metallicity. The advantage of this diagnostics 
is that the lines are closely spaced and easily observed 
in a single spectroscopic setting. However, diagnostics 
using [NII]A6584 are known to depend on the ioniza- 
tion parameter (Kewley & Dopita 2002) and the N/O 
ratio which is not constant with metallicity (e.g. Perez- 
Montero & Contini 2009). These effects may evolve with 
redshift (Kewley et al., in prep). Therefore, an MZ re- 
lation based on this diagnostics may be systematically 
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Fig. 1. — The MZ relation at five epochs ranging to z ~ 2.3. The curves are fits to the data defined by Equation 4. The solid curves 
indicate metallicities determined using the KK04 strong-line method and the dashed curves indicate metallicities converted using the 
formulae of Kcwley & Ellison (2008). Data presented in this figure can be obtained from HJZ upon request. 



biased. Both Erb et al. (2006) and Yabe et al. (2012) 
determine metallicities using the PP04 N2 calibration. 
We convert metallicities determined from the PP04 N2 
diagnostic to the KK04 diagnostic using the conversion 
formulae given in Kewley & Ellison (2008). 

3. RESULTS 
3.1. The MZ Relation 

Figure 1 shows the MZ relations at five epochs. We 
determine the MZ relation for the SDSS, SHELS and 
DEEP2 samples by binning the data. We sort galaxies 
into equally populated bins of stellar mass and plot the 
median stellar mass and metallicity for each bin. The 
MZ relation of Yabe et al. (2012) and Erb et al. (2006) is 
determined from stacked spectra sorted by stellar mass. 
The errors for the z < 1 data are determined from boot- 
strapping. For the z > 1 data the errors are determined 
from the dispersion in the stacked spectra. 

We fit the MZ relation using the function defined by 
Moustakas et al. (2011). The functional form of the MZ 
relation fit is 



12 + log(0/H) = Z a - log 



1 



(4) 



the commonly used quadratic fit (e.g. Tremonti et al. 
2004; Zahid et al. 2011) which turns over at high stellar 
masses. Furthermore, the parameters of the fit reflect 
our physical intuition of chemical evolution and are more 
straightforward to interpret physically (see discussion in 
appendix of Moustakas et al. 2011). In Equation 4, Z a 
is the asymptotic metallicity where the MZ relation flat- 
tens, M is the characteristic mass where the MZ relation 
begins to flatten and 7 is the power law slope of the MZ 
relation for Af* << M Q . The fitted value of Z a is subject 
to uncertainties in the absolute calibration of the metal- 
licity diagnostic, though the relative values are robust 
(see Section 2.3). We do not probe stellar masses where 
M* << M a . Therefore the power law slope of the MZ 
relation at the low mass end, 7, is not well constrained. 
Table 1 lists the fitted parameters. We propagate the 
observational uncertainties to the parameter errors. 



3.2. Scatte 



the MZ Relation 



This function is desirable because it is monotonic unlike 



The scatter in the metallicity distribution as a function 
of stellar mass and redshift provides important additional 
constraints for the chemical evolution of galaxies. In Fig- 
ure 2 we plot the scatter in the MZ relation. We note 
that we have converted the metallicity to linear units for 
clarity. In Figures 2A and 2B we plot the limits contain- 
ing the central 85% and 50% of the galaxy metallicity 
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TABLE 1 
MZ Relation Fit 



Sample Redshift Z \og(M /Mo) 7 Calibration 



SDSS 0.08 9.121 ± 0.002 8.999 ± 0.005 0.85 ± 0.02 KK04 

SHELS 0.29 9.130 ± 0.007 9.304 ± 0.019 0.77 ± 0.05 KK04 

DEEP2 0.78 9.161 ± 0.026 9.661 ± 0.086 0.65 ± 0.07 KK04 

Y12 1.40 9.06 ± 0.36 9.6 ± 0.8 0.7 ± 1.5 PP04 

E06 2.26 9.06 ± 0.27 9.7 ± 0.9 0.6 ± 0.7 PP04 



Note. — The sample and median redshift are given in columns 1 and 2, respectively. The fit parameters from 
Equation 4 are given in columns 3-5. Column 6 indicates the strong line method used for deriving metallicity. 
We convert PP04 metallicities to the KK04 calibration using the formulae from Kewley & Ellison (2008). 




Fig. 2. — Limits containing the central (A) 85% and (B) 50% of the metallicity distribution and the difference between the upper and 
lower limits for the central (C) 85% and (D) 50% of the galaxy metallicity distribution as a function of stellar mass for the z < 1 samples. 
For clarity the metallicity is plotted on a linear scale. 



distribution, respectively, as a function of stellar mass. 
The errors bars are determined from bootstrapping the 
sample distribution. In Figures 2C and 2D we plot the 
scatter in the MZ relation (defined as the difference be- 
tween the upper and lower limits of the 85% and 50% 
contour, respectively) as a function of stellar mass. The 
error bars are the bootstrapped errors on the limits added 
in quadrature. 

4. DISCUSSION 



4.1. Evolution of the MZ Relation 

Figure 1 shows that at a fixed stellar mass (over the 
stellar mass range 10 9 M Q < M* < 10 u Mq) the aver- 
age metallicity of the star-forming galaxy population in- 
creases as a function of cosmic time since z ~ 2.3. How- 
ever, the evolution is mass dependent; thus the shape of 
the MZ relation evolves with time. There is a greater 
metallicity evolution in lower mass galaxies since z ~ 0.8 
as compared to more massive galaxies. The smaller dif- 
ference in the MZ relation for massive galaxies since 
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z ~ 0.8 is due to flattening of the MZ relation at higher 
stellar masses. 

Figures 2A and 2B show that the flattening in the MZ 
relation at late times results from a mass dependent flat- 
tening of the upper envelope of the metallicity distribu- 
tion. The stellar mass where galaxy metallicities saturate 
decreases at late times. The statistically significant evo- 
lution of M D to lower stellar masses at late times (see Ta- 
ble 1) quantifies this effect. The metallicity distribution 
in Figure 2A and 2B also shows that there is an empirical 
upper metallicity limit. Z a in Equation 4 quantifies the 
upper metallicity limit. Within the observational uncer- 
tainties, Z Q is constant to z ~ 2.3. 

Analytical models provide a physical interpretation of 
Z Q and the observed effect of metallicity saturation. In 
a closed-box model of chemical evolution where no gas 
enters or leaves the system, the gas-phase oxygen abun- 
dance can evolve to arbitrarily high values. However, 
galaxies do not evolve as closed boxes. More realistically, 
in an inflow model of chemical evolution where star for- 
mation is fueled by inflowing, pristine gas, the gas-phase 
oxygen abundance saturates at a metallicity equivalent 
to the nucleosynthetic stellar yield 4 (Edmunds 1990). In 
this case, the amount of oxygen produced by massive 
stars is balanced by the amount of oxygen locked-up for- 
ever in low mass stars. When outflows are included in 
analytical models, the upper metallicity limit can -poten- 
tially be reduced to an effective stellar yield (e.g. Ed- 
munds 1990). This yield depends on the composition of 
the outflow, which is not yet well constrained by observa- 
tions. The constancy of Z Q suggests that the metallicity 
of outflows, at a fixed stellar mass, has not evolved sig- 
nificantly since z ~ 2.3. 

Figures 2A and 2B show that the fraction of satu- 
rated galaxies at a fixed stellar mass increases at late 
times. Thus massive galaxies in Figures 2C and 2D have 
smaller scatter in their metallicity distribution at late 
times. The straightforward interpretation of the data is 
that at early times (z ~ 0.8) the scatter in star-forming 
galaxy metallicities is nearly constant as a function of 
stellar mass. As galaxies evolve, metallicity saturation 
occurs first in the most massive galaxies and later at 
lower stellar masses. The stellar mass where the metal- 
licity saturates and the fraction of saturated galaxies at 
a fixed stellar mass evolves with redshift. The data show 
that metallicity saturation leads to a decrease in the scat- 
ter with time (Figures 2C and 2D). Therefore, the ob- 
served scatter in the MZ relation at z ~ 0.8 is likely 
to be a lower limit to the scatter in the MZ relation at 
higher redshifts. 

4.2. Comparison with Previous Observational Studies 

Evolution in the shape of the MZ relation has previ- 
ously been reported (Savaglio et al. 2005; Maiolino et al. 
2008; Zahid et al. 2011). At a fixed stellar mass, there 
is relatively greater enrichment in less massive galaxies 
since z ~ 0.8. Our results provide a robust quantification 
of this evolution. 

Galaxies are observed to form in "downsizing" such 

4 For reference, the saturation metallicity (i.e. gas-phase oxygen 
abundance) in Figure 2 converted to mass density units is ~ 0.013. 
Typical theoretical values of the nucleosynthetic stellar yield of 
oxygen in the same units vary between 0.01 and 0.02. 



that massive galaxies form earlier and more rapidly than 
low mass systems and the star-forming population of 
galaxies is dominated by less massive galaxies at later 
times (Cowie et al. 1996). Flattening of the slope of the 
MZ relation at late times has been interpreted as the 
chemical version of galaxy downsizing (Savaglio et al. 
2005; Maicr et al. 2006; Maiolino et al. 2008; Zahid et al. 
2011). We show that flattening of the MZ relation is 
due to a redshift evolution in the stellar mass where the 
metallicities of galaxies begin to saturate. This is prob- 
ably due to gradual depletion of gas reservoirs at late 
times for galaxies at a fixed stellar mass. In this case, 
the flattening of the slope of the MZ relation can be in- 
terpreted without the need to invoke downsizing. 

Moustakas et al. (2011) study the MZ relat ion to 
z ~ 0.75 using a sample of ~ 3, 000 galaxies drawn 
from the AGN and Galaxy Evolution Survey (AGES, 
Kochanek et al. 2012). In contrast to our results, they 
find that the shape of the MZ relation does not evolve 
in a mass-dependent way for galaxies with stellar masses 
> 10 98 Mq. They base their conclusion on an MZ rela- 
tion which is only derived for the most massive galaxies 
at intermediate redshifts. At z ~ 0.6, the lowest stel- 
lar mass bin has a median stellar mass of 1O 1O 6 M0 as 
compared to 10 9 2 M Q for the DEEP2 galaxy sample at 
z ~ 0.8. Thus the shape of the MZ relation at interme- 
diate redshifts is not well constrained by the data. How- 
ever, a direct comparison of the data suggests that sys- 
tematic differences in sample selection may be the source 
of the discrepancy. The MZ relation measured by Mous- 
takas et al. (2011) at z ~ 0.3 is significantly flatter for 
galaxies with M» > 1O 1O ' 2 M as compared to the rela- 
tion we derive from the SHELS sample at the same red- 
shift. Moustakas et al. (2011) select galaxies with well 
measured [OIII]A5007. Foster et al. (2012) show that a 
S/N cut on the [OIIIJA5007 emission line leads to a flat- 
tening of the MZ relation for massive galaxies because 
metal-rich galaxies have weak oxygen emission and low 
ionization parameter (Zahid et al. 2012a). A systematic 
bias against massive, metal-rich galaxies could lead to 
an underestimate of the relation for massive galaxies and 
may explain the lack of evolution in the shape reported 
by Moustakas et al. (2011). 

5. SUMMARY AND CONCLUSIONS 

We investigate the evolution of the mass-metallicity 
relation and its shape using five samples spanning the 
redshift range of < z < 2.3. We calculate the stellar 
masses and metallicities for all five samples using as con- 
sistent a methodology as is currently possible. Our con- 
clusions arc primarily based on the three z < 1 galaxy 
samples at z — 0.08,0.29 and 0.78 for which we have 
metallicities measured in individual galaxies. We show 
that the MZ relation at z = 1.4 and 2.26 determined from 
stacked spectra are consistent with the main conclusions 
of this study which are: 

• The metallicities of star-forming galaxies at a fixed 
stellar mass decrease at all stellar mass > 10 9 M© 
as a function redshift. 

• The MZ relation since z ~ 0.8 evolves in a mass 
dependent manner such that the shape of the MZ 
relation changes with redshift. 
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Galaxy metallicities saturate. The stellar mass 
where galaxy metallicities saturate and the frac- 
tion of galaxies with saturated metallicities at a 
fixed stellar mass evolves. Thus there is a mass de- 
pendent decrease in the scatter and a flattening of 
the MZ relation at late times. 



• We attribute the flattening of the MZ relation at 
late times to an empirical upper limit in the gas- 
phase oxygen abundance for star-forming galaxies 
which does not evolve significantly. 



We quantify evolution in the shape of the MZ relation 
out to z ~ 0.8. The multiplexing capability of the new 



generation of near- infrared spectrographs will allow us to 
rigorously establish the MZ relation for z > 1 galaxies. 

We thank the referee for useful comments for improv- 
ing the paper, John Moustakas, Vivienne Wilde and Nel- 
son Caldwell for providing useful advice and routines and 
Kiyoto Yabe and Dawn Erb for sharing their data. We 
thank Susan Tokarz, Scan Moran and Warren Brown 
for the careful reduction of the SHELS data and Mike 
Calkins and Perry Berlind for operating Hectospec. HJZ 
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